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Abstract. The structure of the dark matter and the thermodynamical status of the hot gas in galaxy clusters is 
an interesting and widely discussed topic in modern astrophysics. Recently, Rasia et al. (2004) have proposed a 
new dynamical model for the mass density profile of clusters of galaxies as a result of a set of high resolution 
hydrodynamical simulations of structure formation. We investigate the lensing properties of this model evaluating 
the deflection angle, the lensing potential and the amplification of the images. We reserve particular attention to 
the structure and position of the critical curves in order to see whether this model is able to produce radial and 
tangential arcs. To this aim, we also investigate the effect of taking into account the brightest cluster galaxy in 
the lensing potential and the deviations from spherical symmetry mimicked by an external shear. We also analyze 
the implication of the gas density and temperature profiles of the Rasia et al. (2004) model on the properties of 
the X-ray emission and the comptonization parameter that determines the CMBR temperature decrement due 
to the Sunyaev- Zel'dovich effect. 
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1. Introduction 

Being the largest bound structures in the universe, galaxy 
clusters occupy a special position in the hierarchy of cos- 
mic structures in many respects. They can be detected 
at high redshift because of the presence of hundreds of 
galaxies and hot X - ray emitting gas and, therefore, they 
appear to be ideal tools for studying large scale struc- 
ture, testing the theories of structure formation and ex- 
tracting invaluable cosmological information (see, e.g., 
|BorganT fc Guzzo 200TllRosati et al. 20020 . 

To first order, galaxy clusters may be described as 
large dark matter haloes since this component repre- 
sents up to 80% of the cluster mass. Cluster properties 
may thus be investigated by means of numerical simu- 
lations performed in the standard framework of hierar- 
chical CDM structure formation. The impressive growth 
in processing speeds of computers in recent years has 
allowed us to deeply investigate this issue leading to a 
strong debate about this fundamental topic. While there 
is a general consensus that relaxed galaxy clusters ex- 
hibit a density profile that is well described by a dou- 
ble power law with outer asymptotic slope —3, there is 



still an open controversy about the value of the inner 
asymptotic slope (5 with proposed values mainly in the 
range ~ 1.0-1.5 QNavarro et al. 1997IITormen et al. 19971 
IMoore et al. 19981 |Jing fc Suto 20021 IPower et al. 20031 
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INavarro et al. 2004(1 . Un the other hand, a similar contro- 
versy has arisen over the question whether such cusps are 
indeed observed in galaxies (see, e.g., ISimon et al. 2 003 
and references therein). However, on galaxy scale, the ef- 
fect of baryonic collapse and astrophysical feedback pro- 
cesses (such as supernova explosions) may alter signifi- 
cantly the dark halo structure thus complicating the in- 
terpretation of the observations. 

Strongly lensed arcs in galaxy clusters probe the gravi- 
tational potential on scales (r ~ 50—100 kpc) large enough 
to avoid baryonic contamination and are thus an ideal tool 
to investigate this puzzling question. In particular, radial 
arcs probe the slope of the mass profile at their positions, 
while tangential arcs constrain the total mass within their 
radial distance from the cluster centre. Moreover, a mea- 
surement of the velocity dispersion essentially fixes the 
mass divided by the radius even if it is worth stressing that 
this estimate is reliable only for relaxed clusters. These 
observables can then be combined to obtain a powerful 
method to extract information on the cluster structure. 
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It is worth noting that the second most important 
component of a galaxy cluster, namely the gaseous in- 
tracluster medium (ICM), is often neglected in numerical 
simulations. In the usual approach, the ICM distribution 
is determined a posteriori from the dark matter density 
profile imposing the hydrodynamical equilibrium and as- 
suming an isothermal or polytropic equation of state for 
the gas ( |Komatsu fc Seljak 2001) lAscasibar et al. 2003|) . 
However, such an approach is somewhat biased since it 
relies on hypotheses (isothermality and hydrodynamical 
equilibrium) that do not hold in real galaxy clusters. On 
the other hand, it is also possible to determine the radial 
mass profile of both dark matter and ICM directly from 
simulations explicitly taking into account the gas compo- 
nent. This is the approach followed in a recent paper by 
Rasia et al. (2004) . Using an extended set of high resolu- 
tion non radiative hydrodynamic simulations and assum- 
ing spherical symmetry, these authors have first derived 
the phase space density of the dark matter particles and 
then given fitting formulae for the density profile and the 
velocity dispersion thus allowing them to verify the dy- 
namical equilibrium of the system. Turning then to the 
hot gas component, they have derived analytic expressions 
for the density structure, the temperature profile and the 
velocity dispersion of the ICM without imposing any a pri- 
ori hypotheses on the gas dynamical status or its equation 
of state. In particular, Rasia et al. have shown that the 
isothermality hypothesis breaks down at distances from 
the centre larger than ~ 0.2R V , with R v the virial radius 
of the cluster. The Rasia et al. model (hereafter RTM 
model) presents some peculiarities that make it different 
from the other models available in literature. Moreover, 
all the relevant quantities of both the dark matter and 
gas have been derived in a self - consistent way free of any 
bias induced by aprioristic hypotheses on the dynamical 
state of the system. 

In particular, the knowledge of the gas profile allows 
one to resort to a completely different (and complemen- 
tary) observable. With temperature of the order of few 
keV, the ICM gas is dense and hot enough that clusters are 
luminous X - ray sources with the bulk of the X - rays be- 
ing produced as bremsstrahlung radiation f Sarazin "l998). 
Electrons in the ICM are not only scattered by ions, but 
may themselves Compton scatter photons of the cosmic 
microwave background radiation (CMBR) giving rise to 
the Sunyaev- Zel'dovich (SZ) effect (see IBirkinshaw 1999 
for a comprehensive review). The temperature decrement 
due to the SZ effect is able to provide information on 
the cluster structure, on the motions of galaxy clusters 
relative to the Hubble flow and on the Hubble flow it- 
self and the cosmological constants that characterize it 
l|Reese et al. 20021 ISereno 2003|). 

In order to investigate if the RTM model is a viable 
one, a direct comparison with the main cluster obscrvables 
(both from lensing and SZ effect) is needed. As a first step, 
one has to study the lensing properties of the RTM model 
and to calculate the SZ effect taking care of the peculiar 
temperature profile. This is the aim of the present paper. 



In Sect. 2 we evaluate the deflection angle and the lens- 
ing potential of the spherically symmetric RTM model. 
Radial and tangential arcs form near the position of the 
critical curves. Therefore, Sect. 3 is devoted to a detailed 
investigation of the critical curves structure of the model 
with a particular emphasis on how these properties de- 
pend on the model parameters. The effect of taking into 
account the contribution of the brightest cluster galaxy 
to the lensing potential is investigated in Sect. 4, while 
the impact of deviations from spherical symmetry or tidal 
perturbations from nearby clusters are mimicked by an 
external shear and discussed in Sect. 5. Having been pro- 
posed recently, the RTM model has to be compared with 
the previous proposals. In particular, in Sect. 6, we com- 
pare some of its lensing properties with those of the NFW 
model investigating possible systematic errors in the virial 
mass estimate. The computation of the SZ effect due to 
the distribution of the ICM in the RTM model is presented 
in Sect. 7, while the results are compared to the prediction 
of both the (3 model and the NFW model in Sect. 8 where 
we also discuss the detectability of RTM clusters in SZ 
survey. The details of the numerical simulations on which 
the RTM model is based may induce systematic errors on 
the main results. Some qualitative comments on this topic 
are presented in Sect. 9. We summarize and conclude in 
Sect. 10. 

2. Deflection angle and lensing potential 

Let us adopt a rectangular coordinate system (x, y, z) 
with origin in the cluster centre and let (r, 9, 4>) be the 
usual spherical coordinates. The mass density profile of 
the RTM model is IjRasia et al. 2004|l : 



p{r) = PoPb 
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with pi, = ^luPcrit the present day mean matter density 
of the universe and : 



Pa 



(1 - fb)A v 



6[(l + 2x p )/(l + x p y/ 2 -2x 1 p /2 } 
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where A„ is the virial overdensity specified by the cos- 
mological model and the term is the average baryonic 
fraction used to properly weight the dark matter compo- 
nent in the cluster. Following Rasia et al. (2004), we set 
fb = 0.097 as obtained by averaging over their sample of 
simulated clusters. 

The RTM model is fully characterized by two param- 
eters, namely the dimensionless scale radius x p (or the 
concentration crtm = l/^p) an d the virial radius R v . 
However, it is more convenient to express R v in terms of 
the virial mass (i.e. the total cluster mass) M v using the 
following relation : 



R v 



( 3M„ 
\AirA v p b 



1/3 



(3) 
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We will assume that the model is spherically symmetric 
so that all the lensing quantities will depend only on the 
projected radius R = (x 2 + y 2 ) 1 ! 2 . This is the same ap- 
proximation used in Rasia et al. (2004) to obtain the den- 
sity profile in Eq.QJ. While useful in the computations, 
this approximation is not a serious limitation to our anal- 
ysis since the results for the circular case may be immedi- 
ately generalized to flattened models by means of numer- 
ical integration l|Schramm 19901 iKeeton 20"oT|l . Moreover, 
we will investigate later the impact of deviations from cir- 
cular symmetry by adding a shear term to the lensing 
potential. However, we stress that the results for the cir- 
cularly symmetric models allow us to obtain a picture of 
the main properties of the RTM model as a lens. 

As a first step to investigate the lensing properties of 
the RTM model, we have to evaluate the corresponding 
surface mass density. Starting from the definition : 
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Fig. 1. The scaled deflection angle a norm = a/a v vs the 
dimensionless radius £ for three values of the dimensionless 
scale radius x p of the RTM model, i.e. x p = 0.1 (solid), 
x p = 0.2 (short dashed), x p = 0.3 (long dashed). 



p(x,y,z)dz 



and using a convenient transformation to spherical coor- 
dinates, we get : 



£(£) = 2 p p b R v / 
Jo 
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with £ = R/R v and £„ the surface density at the virial 
radius given by : 



E(£ = l) 



p Q p b R v S(l,x p ) 



(5) 



and wc have defined the function : 
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In the previous equation, r(£) is the actual T function and 
pF g [{ai, ...,a p }, {h, . . . , b q }, () is the generalized hyper- 
geometric function 1 ( |Gradshteyn fc Ryzhik 1980| ). 

Having obtained the surface mass density, it is now 
straightforward to compute the deflection angle. Because 
of the circular symmetry in the lens plane, the de- 
flection angle is purely radial and its amplitude is 
(jSchneider et al. 1992|l : 



a= R 



2 ^S^W 



(7) 



1 We use the same notation for the generalized hypergeomet- 
ric function as in the Mathematica package. 



with E cr it = c 2 D s /4irGDiDi s the critical density for lens- 
ing and D s , Di, Di s are the angular diameter distances 
between observer and source, observer and lens and lens 
and source, respectively. Inserting Eq.QJ into Eq.Q, we 
get: 



"(0 = a v x 



(8) 



F(l,x p ) 

with a v the deflection angle at the virial radius given by 
AY, v R v T{l,Xp) 
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and we have introduced the function : 
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It is interesting to observe that only the parameter x p of 
the RTM model determines the behaviour with the dimen- 
sionless radius £ of the scaled deflection angle a/a v , while 
the virial radius R v (and thus the total mass M v ) enters 
only as a scaling factor. As an example, Fig. ^ shows a/a v 
for three values of x p . This plot may be qualitatively ex- 
plained as follows. Lower values of x p correspond to higher 
values of the concentration crtm and hence to more mass 
within a fixed radius. Since, for a given position in the lens 
plane, the deflection angle scales with the projected mass 
within £, it turns out that a/a v is higher for lower values 
of Fig. n shows. 

This qualitative discussion also explains the behaviour 
of a v with Xp that is shown in Fig. [21 A comment is in 
order here to understand how this plot has been obtained. 
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Fig. 2. The deflection angle at the virial radius a v vs the 
dimensionless scale radius x p of the RTM model. See the 
text for the values of the other parameters. 



Eq.Q shows that a v depends on the cluster parame- 
ters (M„, x p ), the lens and source redshift (zi, z s ) and the 
background cosmological model. We adopt a flat ACDM 
model with (» M , ^A, h) = (0.3, 0.7, 0.72) giving A v = 324 
l(Eke et al. 199611 . We set {z h z s ) = (0.313,1.502) as for 
the real cluster lens MS2137-23 (Sand et al. 2002, 2004). 
Unless otherwise stated, the same values for the cosmo- 
logical parameters and the lens and source redshift will be 
used throughout the paper. To obtain the plot in Fig.[2I 
we have fixed M v = 7.5 x 10 14 M Q , but the results for 
other values of M v may be easily scaled observing that 
a v oc T, V R V oc R 2 oc M^ 3 . 

Let us now derive the lensing potential ip(R) for the 
RTM model. To this aim, one should solve the two dimen- 
sional Poisson equation ((Schneider et al. 1992(1 : 



V 2 ip 



2k 



(11) 



with k = S/E cr i t the convergence. However, because of 
the circular symmetry in the lens plane, it is also : 

a(R) = ^ -> MR) = [ a(R)dR . (12) 
dR J 

Inserting Eq. (JHJ into Eq. 1(12(1 and integrating, we find : 



with : 
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Finally, let us consider the lens equations. Adopting po- 
lar coordinates (R, in the lens plane with $ measured 
counterclockwise from North, the time delay of a light ray 
deflected by the gravitational field of the cluster lens is : 



At = h- 1 



Two x 



-R 2 - RR S cos (tf - tf s ) + -R 2 S - ip(R, i?) 



(16) 



where (i?,i9) is the image position, (_R s ,"i? s ) the unknown 
source position and we have defined : 



Two = 



DiD s \ (1 + jbi) 
Dis 



(17) 



According to the Fermat principle, the images lie at the 
minima of At, so that the lens equations may be simply 
obtained by minimizing At. Inserting Eq.O into Ea. lfTfi)) 
and differentiating, we get : 



£-&cob (0-4.) 



e s sin(^-i? s ) = 



R v T(l 7 x p ) 



(18) 



(19) 



with ^ s = R s /R v . En. 1)19(1 has two solutions. The first is 
£ s = 0, i.e. lens, source and observer are perfectly aligned. 
The only image is the Einstein ring that we will discuss 
in much detail in the next section. The second solution 
is obtained for sin t9 — i9 s = <^=> 1) = d s + rmr with 
m = 0, 1. In this case, we get two images 2 symmetrically 
placed with respect to the lens centre. The radial coordi- 
nate £1 of the first image (i.e., the one with d = $ s ) is 
obtained by solving Eq. 1(18(1 with cos($ — i? s ) = 1, while 
the second has a distance £2 from the lens centre obtained 
by solving Eo. 1(18(1 with cos (1? — $ s ) = —1. These equa- 
tions may be solved numerically provided that the clus- 
ter parameters have been fixed. Up to now, there is only 
one multiple image system in which the lens is a cluster 
galaxy, namely the recently discovered SDSS J1004+4112 
l|Inada et al. 20031 |Oguri et al. 2004D , while for all the 
other multiply imaged quasar the lens is a galaxy (see, e.g., 
the CASTLES web page. lKochanek et al. 20031. Since the 
RTM model has not been tested on galactic scale (given 
the mass range probed by the simulations employed by 
Rasia et al. 2003), we prefer to not discuss further the 
formation of multiple images. 

3. Critical curves 

The most spectacular effect of lensing by a galaxy cluster 
is the formation of giant arcs (see, e.g,. IKneib et al.~1 996 
for the textbook example of A2218). These are very lumi- 
nous and highly distorted images of a source galaxy whose 



-2yj2i:Xp ln£ 



Note that the theorem of the odd number of images 
ISchneider et al.~ 992 Strauma nn 19981 does not apply here 
(15) because the mass distribution is singular in the origin. 
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Fig. 3. Contours of equal Einstein radius Re in the 
(log M v ,Xp) plane. Re ranges from 5 (the uppermost 
curve) to 35arcsec (the lowermost one) in steps of 
10 arcsec. 
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Fig. 4. Contours of equal radial critical curve distance 
Rrad from the cluster centre in the (logM v ,x p ) plane. 
Rrad ranges from 2.5 (the uppermost curve) to 5.5 arcsec 
(the lowermost one) in steps of 1 arcsec. 



position is near one of the critical curves of the lensing po- 
tential. The position of the arcs in a given lensing system 
allows to strongly constrain the cluster mass distribution 
and can also be used to determine the cosmological param- 
eters IjSereno 2002fl . Moreover, arc statistics is a promising 
and efficient tool to discriminate among different cosmo- 
logical models and theories of structure formation. It is 
thus quite interesting to investigate the number (and the 
type) of arcs the RTM model may form. 

To this aim, let us first remember the expression for 
the magnification /x of a source due to the lensing effect. 
It is IjSchneirler et al. 1992|l : 



1 



1 



1 



M = 



detA (1 - ip xx ){l - tpyy) - i>l v X r Xt 



(20) 



where A is the amplification matrix (that is the jacobian 
matrix of the lens mapping) and (A r , At) for a circularly 
symmetric model are given as : 



A '' 1 dR ' 



A, 



1 



a 
R 



(21) 
(22) 



Inserting Eq.JSJl into Eqs.J2U, we get the corre- 

sponding quantities for the RTM model that we do not 
explicitly report here for sake of shortness. The critical 
curves are the loci where det A = 0. This condition is sat- 
isfied by imposing A r = or At = 0. The second equation 
implicitly defines the tangential critical curves. It is easy 
to show that : 



A t = 



R v T(l,x p ) 



(23) 



with (,e = Re/Rv and Re the Einstein radius. Ea. (|23[l 
may be solved numerically for fixed values of the clus- 
ter model parameters. Fig.|31 shows the contours of equal 
Re in the (\ogM v ,x p ) plane 3 with higher values of Re 
corresponding to lower curves in the plot. We have only 
considered clusters with \ogM v £ (14.5, 15.5) because this 
is (approximately) the mass range probed in Rasia et al. 
(2004), while we will (usually but not always) consider 
x p e (0.01/7.13,2/7.13) since Rasia et al. (2004) states 
that the average value of (crtm) — (l/^p) — 7.13 over 
their cluster sample. 

Fig-EH shows that highly concentrated clusters (i.e., 
with small values of x p ) give rise to tangential critical 
curves that are more distant from the cluster centre. Thus, 
one could qualitatively conclude that only RTM models 
with lower values of x v are able to produce tangential 
arcs. However, one should also take into account that, as 
expected, for a fixed x p the Einstein radius increases with 
M v . Therefore, RTM models with high values of x p could 
still produce tangential arcs with large radii provided that 
the mass is large enough and that it is possible to extrap- 
olate the RTM model outside the mass range probed by 
the simulations. 

As said above, radial critical curves are implicitly de- 
fined by the condition A r = 0. Inserting Eq.Q and JSJ 
into Eq.(j2U, we get R ra d, the radial critical curve dis- 
tance from the cluster centre. To investigate how R ra d 
depends on the model parameters, we plot in Fig.01 the 
contours of equal R ra d in the (log M v ,x p ) plane. The be- 
haviour of R rac i with the log M v is qualitatively similar to 



Unless otherwise stated, log x is the logarithm base 10. 
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that of the Einstein radius Re, but the opposite holds for 
the dependence on x p . For a given value of logM„, R ra d 
increases with x p , while Re decreases. As a result, highly 
concentrated RTM models (i.e. models with small x p and 
henceforth high crtm) give rise to tangential arcs situated 
at large distances from the cluster centre, but the radial 
arc lies in the very inner regions of the cluster. 

As a simple application of these results, we consider 
the case of the real cluster lens MS2137-23 in which both 
a tangential and a radial arc have been observed with 
Re = 15.35 arcsec and R ra d = 4.5 arcsec (Sand et al. 2002, 
2004). An easy way to find out the values of the RTM 
model parameters able to fit the arc positions in MS2137- 
23 is to draw the contour levels in the plane (log M v , x p ) 
for Re and R ra d equal to the values quoted above and 
look for the intersection point between these two curves. 
It turns out that the best fit parameters are (logM„, x p ) = 
(14.7, 0.019), i.e. M v = 5.0 x 10 14 M and c RTM ^ 53. As 
a further test, we also consider the case of the real cluster 
lens RXJ1133 at redshift z L = 0.394 IjSand et al. 2004|) . 
Both a radial and a tangential arc are observed with z$ = 
1.544 and (R ra d, Re) = (3.2, 10.9) arcsec that may be ob- 
tained by describing the cluster with a RTM model with 
best fit parameters (\ogM v ,x p ) = (14.5,0.018), i.e. M v = 
3.2 x 10 14 M Q and crtm — 56. The values of the concen- 
tration are quite high if compared to (crtm) — 7.13 found 
by Rasia et al. (2004) over their sample of simulated clus- 
ters. However, this could not be considered an evidence 
against the RTM model. Actually, we have only considered 
the spherical case, while it is well known that also a small 
cluster cllipticity changes significantly the position of the 
critical curves (see, e.g., |Bartelmann fc Meneghetti 20 03). 
Moreover, one should also take into account the impact on 
the critical curves of the galaxy lying at the centre of the 
cluster gravitational potential and of other eventual sub- 
structures. That is why we do not speculate further on the 
high crtm values needed to reproduce the arcs positions 
in MS2137-23 and RXJ1133, while a detailed comparison 
with observations will be presented elsewhere. 

4. Adding a bright cluster galaxy 

Some recent studies have highlighted the importance of 
considering the brightest cluster galaxy (hereafter BCG) 
when investigating the lensing properties of a cluster 
< |Meneghetti et al. 2003|lSand et al. 2004|) . It is thus inter- 
esting to study how the critical curves of the RTM model 
are affected by the addition of a BCG. To this aim, we 
place the galaxy exactly at the centre of the cluster and 
model it using the Hernquist profile whose mass density 
is flHernquist 1990| ) : 



p{r) 



Ps 



r/r s (1 + r/r s ) 3 



(24) 



with p s a characteristic density and r s a scale radius. 
The Hernquist profile has the notable property that 
its projected density well approximates the i? 1 / 4 law 
Ijde Vaucouleurs 1948|l provided that the effective radius 
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Fig. 5. Zero level curves of the denominator in Ea. (|28|l 

in the plane (R,\ogM v ) for a RTM cluster with x p = 
0.5/7.13. The solid line refers to the case with no BCG, 
while the dashed one has been obtained adding a BCG 
with (M h ,R e ) = (5 x 10 12 M Q , 24.80 kpc). For a given 
value of M v , the intersections of the horizontal line 
logM„ = const with the curves plotted are the values 
of R ra d and Re respectively. 

is related to the scale radius of the Hernquist model by 
the relation: R e ~ 1.81r s . The model is fully character- 
ized by two parameters that we choose to be R e and the 
total mass Mh given by : 

M h = 2-kt\ Ps . (25) 

Assuming spherical symmetry, the deflection angle of the 
galaxy is IjKeeton 2001fl : 



a H {<j) = 2k s t, 



a 2 - 1 



(26) 



with a = R/r s and k s — p s r s /Y, cr i t and we have defined : 

(27) 



^=|= tan- 1 \/^2^T (a > 1) 
H((t) = <j 7 Ji_tanh- 1 VT — ^ ((T < !) 

1 ° (a = l) 



Since the standard theory of lensing is developed in the 
weak field limit, the total deflection angle is simply the 
sum of the contributions from the cluster and the BCG. 
Hence, the total magnification may be written as : 



1 



/' 



{\h + \?™ - l)(Af + \f™ - 1) 



(28) 



with A r and At given by Eas. (|21|) and 122|) respectively and 
quantities with the superscript "H" ( "RTM") refers to the 
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Fig. 6. Same as Fig.[5]but now the value of M v is fixed as 
1.125 x 10 15 M© and x p is changing. Solid line refers to 
the case with no BCG, while dashed one to the case with 
a BCG with parameters fixed as before. 



Hcrnquist (RTM) model. The radial and tangential arcs 
radii R ra d and Re are defined as those radii vanishing the 
first and the second term respectively of the denominator 
in Eq.{2SJ. 

To study where the critical curves form, we plot in 
Fig. |S] the loci in the (R, log M v ) plane where the total 
magnification formally diverges having arbitrarily fixed 
x p = 0.5/7. 13 and R e = 24 .80 kpc as for the BCG in 
MS2137-23 l|Sand et al. 20040 . The solid line refers to the 
case with no BCG, while the dashed one shows how the 
curves are modified by the addition of a BCG with total 
mass 4 Mh = 5 x 10 12 M©. The number of critical curves is 
still two, but their position is affected by the presence of 
BCG with the distance between them increased with re- 
spect to the case with no BCG. In particular, while R ra d is 
slightly smaller or higher depending on the ratio between 
the mass of galaxy and that of the cluster, the Einstein 
radius Re significantly increases. This is expected since 
Re is proportional to the total mass within the tangen- 
tial critical curve so that, adding the BCG mass, Re, gets 
obviously higher. 

Fig.Elis similar to Fig.0 but now the cluster mass is 
set as M v = 1.125 x 10 15 M© and we let x p changing, while 
the BCG parameters are fixed as before. The increase of 
Re is still visible, but what is most important to note 
is the possibility to have radial and critical curves with 



4 Note that the value chosen for the mass of the BCG, 
even if high, is not unrealistic. For instance, the esti- 
mated total mass of the Milky Way is l-9±f;? X 10 12 M 
(Wilkinson & Evans 1999) which is in the mass range we have 
adopted for the BCG. 
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Fig. 7. Constraints in the (log M v , x p ) plane imposing 
Rrad = 4.5 arcsec (left) or Re = 15.35 arcsec (right) for 
the RTM model with (dashed line) and without (solid line) 
taking into account the presence of a BCG. The galaxy pa- 
rameters are set as {Mh, R e ) = (5.0 x 10 12 M Q , 24.80 kpc). 

higher values of x p . From Fig.|^l one sees that, when the 
BCG is absent, RTM models with x p > 0.18 are unable to 
produce radial arcs (i.e. it is R ra d = 0), while the radial 
critical curve appears when the BCG is taken into account 
in the total lensing potential even for less concentrated 
(i.e. with larger x p ) clusters. 

Finally, we investigate qualitatively how the con- 
straints on the RTM model parameters are changed by 
the presence of a BCG. To this aim, we plot in Fig.[3 
the constraints imposed by the presence of a radial arc at 
Rrad = 4.5 arcsec or a tangential arc at Re = 15.35 arcsec 
as observed for the real cluster lens MS2137-23 (Sand et 
al. 2002, 2004). The results in the right panel are easy to 
explain qualitatively. Adding a BCG pushes up the curve 
in the (log M v , x p ) plane so that a tangential arc at a given 
distance may be produced by less concentrated and less 
massive clusters with respect to the case with no BCG. As 
yet noted, Re is proportional to the total projected mass 
within Re itself. Since now the BCG provides part of this 
mass, less mass has to be contributed by the cluster and 
thus less massive and concentrated models are needed to 
obtain a given value of Re- Note, however, that the de- 
viations from the case with no BCG are quite small as 
expected given the high mass ratio between the cluster 
and the galaxy. On the other hand, the left panel shows 
that adding a BCG requires more concentrated and mas- 
sive clusters to produce a radial arc at a given R ra d with 
respect to the case with no BCG. It is also worth noting 
that the constraints from the position of the radial arc are 
more sensitive to the presence (or absence) of the BCG 
(see, e.g., the distance between the solid and dashed line 
in the left panel compared to the same in the right one). 
This is expected since the radial critical curve is inner- 
most and thus probes a range that is more sensitive to 
the inner structure of the cluster where the BCG plays a 
more significant role. 



5. The impact of the shear 

Up to now, we have investigated the lensing properties 
of the RTM model (with and without a central BCG) 
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Fig. 8. Critical curves (solid) and caustics (dashed) for 
the RTM model with the addition of an external shear. 
Upper panels are for 7 = 0.0 (left) and 7 = 0.05 (right), 
while lower panels refer to 7 — 0.10 (left) and 7 = 0.15 
(right). The shear position angle is set at 9 1 = 0, while the 
cluster parameters are (M v , x p ) = (1.125 x 10 15 M©, 0.03). 

assuming spherical symmetry of the mass distribution. 
However, real clusters are moderately elliptical and it is 
well known that even small ellipticities may alter sig- 
nificantly the lensing properties of a given model. In 
particular, taking into account deviations from spheri- 
cal symmetry is very important when trying to extract 
constraints on the model parameters from the position 
of the lensed arcs in real systems (as clearly demon- 
strated, for instance, in |Bartelmann fc Meneghetti 2 003 
Dalai & Kceton 2003). On the other hand, even when as- 
suming spherical symmetry, it is important to take into 
account also the effect of substructures in the cluster mass 
distribution and possible tidal deformations due to nearby 
clusters. Finally, the large scale structure as a whole could 
also have a not negligible effect i|Keeton et al. 19 97). 

To the lowest order, all these effects may be mimicked 
by adding an external shear to the lensing potential which 
is now written as : 

i>(R,§)=Tj; RTM (R)-^jR 2 co S 2(#-# y ) (29) 

with iPrtm given by Eq. l|13|) and (7,$ 7 ) the shear 
strength and position angle. Without loss of generality 
we assume that the shear is oriented along the major axis 
so that it is fl 1 — 0. Fig. [H] shows the critical curves and 
the caustics for the lensing potential in Ea. (|29|) . To be 
quantitative, we also report in Table 1 the quantity : 
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0.05 
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13 


13 


-12 


0.10 


-22 


29 


28 


-22 


0.15 


-31 


46 
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-31 



Table 1. The impact of the shear on the length along the 
major an minor axes of the tangential and radial critical 
curves. The upper half of the table refers to RTM model 
with x p — 1/(crtm), while for the lower half it is x p = 
0.5(crtai)- In both cases, the virial mass is set to M v — 
1.125 x 10 15 M R . 



with £i the length of the radial or tangential critical curve 
along the i- axis (with i = x,y). 

In the case with no shear (upper left panel), the two 
critical curves are spherical, while the tangential caustic is 
the origin and the radial one is a circle. The effect of the 
shear is to deform the critical curves into ellipses, while 
the inner caustic (corresponding to the tangential criti- 
cal curve) takes a diamond shape and the external one 
becomes elliptical. In particular, the radial critical curve 
is more and more elongated along the major axis 5 as the 
shear strength increases. Quantitatively, this could be seen 
in Table 1 where A£ a is positive and increases with 7, 
while A£ y is negative and different from A£ x . This is what 
corresponds graphically to an ellipse more and more elon- 
gated along the major axis as 7 increases. On the other 
hand, the tangential critical curve may also deviate from 
the elliptical symmetry taking a dumbell shape as 7 gets 
higher. As a consequence, the higher is 7, the higher is the 
the radial arc distance (measured along the major axis), 
R ra d, from the cluster centre for fixed values of the RTM 
model parameters. A similar result holds also for Re- It is 
worth noting that these effects are qualitatively the same 
whatever is the value of x p , but are more pronounced (i.e. 
Mi is larger) for lower values of x p (i.e. higher concentra- 
tions) as can be seen from Table 1 . 

Note that these results are in qualitative agreement 
with the approximate analytical treatment presented in 
Bartelmann & Meneghetti (2004). Actually, these authors 
used a different approach deforming the lens model so that 
the isocountour lines of the lensing potential are ellipses 
with cllipticity e. Elliptical deformation of the lensing po- 
tential leads to dumbell shaped surface mass distribution 
for values of e > 0.2. Even if clusters are highly struc- 
tured, similar mass models are quite unrealistic so that we 
have preferred not to follow this approach. On the other 
hand, it is possible to show that an elliptical potential 
ip{x 2 + y 2 /q 2 ) with an on axis shear and axial ratio q pro- 
duces the same image configuration as a pure elliptical 



5 This is a consequence of having fixed ■d 1 — 0. In general, 
the ellipse is elongated along the direction individuated by $ 7 . 
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Fig. 9. lZ a vs £ for the average values of the concentration 
parameters of the NFW and RTM model. 



potential with axis ratio q' = q\J (1 — 7)/(l + 7) without 
shear (|Witt 19960 . In our case, this means that using the 
lensing potential given by Eq. H29|) is equivalent to deform 
the RTM model such that the lensing isopotcntial con- 
tours are ellipses with axis ratio q' = y/ (1 — j)/(l + 7). 
This shows the complete equivalence among our approach 
and that of Bartelmann & Meneghetti (2004). 

6. Comparison with the NFW model 

It is interesting to compare the lensing properties of the 
RTM model with those of the model proposed by Navarro, 
Frenck & White (1997, hereafter NFW) and mostly used 
in literature. Using the same normalization as in Rasia et 
al. (2004), the density profile of the NFW model is : 



ft 



P0,NFWPb 



(x/x s )(l + x/x s ) 2 



(30) 



with x = r/R v , x s = 1/cnfw, cnfw the concentration 
of the NFW model and : 



P0,NFW 



(1 - h)Av 



3 [In (1 + c NFW ) - cnfw/ (1 + c NF w)} 



(31) 



The deflection angle for the NFW model may be conve- 
niently written as (Bartelmann 1996, lKeeton 2 0011 : 



„NFW 



In (c NFW £/2) + TL{c N fw€) 
In (cnfw/2) + 'H{cnfw) 



(32) 



having defined : 



NFW 



4(l-A)A^ x ^ x 
3S cr a 

In (cnfw/2) + H{cnfw) 
In (1 + cnfw) — cnfw 7 (1 + cnfw) 



with R v expressed in arcsec 

Using Eqs.©, ©, © and J33J| and the definition of 
virial radius, we get the following expression for the ratio 
between the deflection angle of the NFW and RTM model : 
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Fig. 10. lZ a vs £ setting the NFW concentration to the 
value predicted by Eq.^SJ and crtm — 7.13 for three 
values of the virial mass, namely M v — 5 x 10 14 (short 
dashed), M v = 7.5 x 10 14 (solid), M v = 10 15 (long 
dashed) . 



CtNFW 
CiRTM 



A78 / M™ \ 1/3 

e v Mv ) 



(1 + 2x p )/y/l + x p — 2^Jx~ p 



ln(l + c NFW ) 
In (c NFW £/2) - 



- cnfw /(I + cnfw 

H(CNFwi) 



(34) 



It is worth noting that the ratio does not depend on the 
redshift of lens and source as it is expected since it is re- 
lated to the different density profiles of the two models 
which is, of course, the same at all redshifts. Moreover, 
to better compare the lensing properties of the two mod- 
els, it is meaningful to assume that the virial mass is the 
same so that, at a given radius £, lZ a only depends on the 
concentrations crtm and cnfw of the two models. Fig.|^l 
shows 1Z a (0 setting crtm = 7.13 and cnfw = 6.8 as de- 
termined by Rasia et al. (2004) averaging over their sim- 
ulated clusters sample. The deflection angle of the NFW 
model turns out to be higher than that of the RTM model 
until £ < 0.6, while lZ a < 1 for light rays impacting in 
the outer region of the halo. This result is simply re- 
lated to the different mass profile of the two models with 
Mnfw {O/Mrtm (£) being larger than 1 in the inner re- 
gions for this choice of (cnfw, crtm)- 

According to some authors, the NFW model is a one 
parameter model since it is possible to relate the concen- 
tration cnfw to the virial mass even if this relation has 
a quite large scatter. Following Bullock et al. (2001), we 
adopt : 



(33) cnfw — 15 — 3.3 log 



W 2 h- 1 M 







(35) 



and plot, in Fig. llQI lZ a for three different values of the 
virial mass. The qualitative behaviour is the same, but it 
is more pronounced for lower mass models corresponding, 
according to Ea. l|35|) . to lower concentrations. 
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Fig. 11. Contour plot in the (logM„, cnfw) plane for the 
Einstein radius of the NFW model. Left panel: contours 
corresponding to values of Re for the RTM model with 
Xp = l/(c RTM ) and M v from 7.5 to 17.5 x 10 14 M (from 
left to right) in steps of 2.5 x 10 14 M Q . Right panel: con- 
tours corresponding to values of Re for the RTM model 
with M v = 10 15 and x p from 0.06 to 0.14 (from right to 
left) in steps of 0.02. The dashed line refers to Ea. l|35|) . 



The NFW model is used in most of the studies of the 
structure of dark matter haloes in galaxy clusters. It is 
therefore interesting to investigate what is the error in- 
duced by using the NFW model to describe a cluster 
that is actually better described by the RTM model. As a 
straightforward example, we consider the estimate of the 
virial mass from the size of the Einstein radius. To this 
aim, let us look at Fig. ^2 where we plot the contour level 
curves in the (logM„, cnfw) plane for the Einstein radius 
of the NFW model corresponding to values of Re evalu- 
ated for RTM models. Let us consider, for instance, the 
RTM model with (M v ,x p ) = (1.5 x 10 15 M , 1/7.13) for 
which one obtains the fourth line (from left) in the left 
panel of Fig.rH] If one assumes that Ea. i|35)l holds, then 
one should estimate the NFW model parameters from the 
intersection point of the dashed curve with the fourth line 
thus grossly underestimating the mass. Actually, even if 
one does not use Ea. l|3"5|l . M v turns out to be underesti- 
mated for the values of cmfw in the plot. To get the cor- 
rect value of M v , one should select a value of crtm that 
is unrealistically low for a galaxy cluster. We thus con- 
clude that using the NFW model to study a cluster that 
is intrinsically described by the RTM model leads to un- 
derestimate the virial mass by an amount that depends on 
the concentration of both the NFW and the RTM model. 

In principle, one should compare the lensing properties 
of the NFW and of the RTM models by also including the 
effects of the brightest cluster galaxy and of the shear. 
However, this should increase the number of parameters 
to eight : two for the NFW model, two for the RTM model, 
the BCG mass and scale radius and the shear strength 
and orientation. It is likely that some degeneracies could 
occur among parameters rendering a comparison between 
the NFW and the RTM model meaningless in this case so 
that we prefer to not perform this test. 



7. The thermal Sunyaev- Zel'dovich effect 

At the beginning of Seventies, Sunyaev and Zel'dovich 
(1970, 1972) suggested that the cosmic microwave back- 
ground radiation (CMBR) can be scattered by the trapped 
hot intracluster electrons giving rise to a measurable dis- 
tortion of its spectrum. This (inverse) Compton scattering 
(now referred to as Sunyaev - Zel 'dovich effect, hereafter 
SZE) has been recognized in the last two decades as an 
important tool for cosmological and astrophysical stud- 
ies (Birkinshaw 1999). More recently, the SZE has been 
used to investigate several physical properties of the gas 
with much attention devoted to the geometry of its den- 
sity profile as well as its thermodynamical status. In par- 
ticular, it has been shown that the canonical hypotheses 
of spherical symmetry and isothermal temperature profile 
may induce errors (of order up to 30%) on the estimated 
values of different parameters such as the Hubble constant 
dRoettiger et al. 1997| IPiffaretti et al. 2003|l . 

The physics of the SZE is quite simple to understand. 
A gas of electrons in hydrostatic equilibrium within the 
gravitational potential of a cluster will have a tempera- 
ture T e ~ GMm p /(2k B Reff) ~ few keV, with M and 
Ref f typical values of the total mass and size of a clus- 
ter. At this temperature, the thermal emission in X-ray 
is composed of thermal bremsstrahlung and line radiation 
processes. Electrons in the intracluster gas are not only 
scattered by ions, but can themselves scatter photon of 
the CMBR giving rise, in average, to a slight change in 
the photon energy. Because of this inverse Thomson scat- 
tering, an overall change in brightness of the CMBR is 
observed. As a consequence, the SZE is localized and vis- 
ible in and towards clusters of galaxies having an X-ray 
emission strong enough to be detectable. 

In the non relativistic limit, the scattering process can 
be described by the Kompaneets equation 



dn 
dy 



P. 



1__9_ 

■ 2 dp e 



On 
dp, 



(36) 



which describes the change in the occupation number of 
photons n{y). In Eq. (|36|l . we have introduced the dimen- 



sionless variable 6 



Pe 



hv e /kBT e , while y is the so called 



Comptonization parameter defined as : 



y 



n e o~T<ll 



(37) 



being <jt the Thompson scattering cross section and n e 
the electrons number density and the integral is performed 
along the line of sight. Assuming that the photons distri- 
bution after the scattering is close to the equilibrium one, 
we have that : 



dn 
dp, 



> n, n 



(38) 



Here, we denote with the underscript "e" quantities refer- 
ring to electrons, while the same quantities without underscript 
refers to photons. 
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so that Ea. l3(jfl simplifies to 



7.1. The gas profile and the structure integral 



dn 
dy 



p 2 ^ dp 



(39) 



having also replaced p e with p because of the homogeneity. 
By solving this equation in the quasi equilibrium hypoth- 
esis, it is possible to obtain both the variation An in the 
occupation number with respect to the equilibrium value 
no and the corresponding shift in temperature AT. It is : 



An y p e v 
n 



eP-1 



pcoth 



-4 



AT 



pcoth ( -J - 4 =y g{p) , 



(40) 
(41) 



where g(p) is the SZE frequency spectrum and we have 
considered that T e (~ 10 7 K) is much higher than the 
CMBR temperature T ~ 2.7K. In the limit of low fre- 
quencies, we get the useful approximated expression : 



AT 



(42) 



Eq. 1)42(1 allows to evaluate the shift in temperature due to 
the SZE provided that the gas number density n e (r) and 
the temperature profiles T e (r) are given. For the RTM 
model, it is IjRasia et al. 2004|) : 



ne{s) = P9fiPb {s + Xp) -2. 5 



Te(s) 



/J, m p 

rp rp 0.016 

1q1 v S 



0.13 



(43) 
(44) 



where s = r/R v , fj, is the mean molecular weight, m p the 
proton mass, T v the virial temperature and p Si o & normal- 
ization density given by : 



Pg,o 



fb&v 



2 + I0x p + (40/3)4 + ( 16 /3) x l 16 1/2 



(45) 



Finally, in Eas. l43|l and (^4 If) . x p and To are fitting param- 
eters to be determined on a cluster by cluster basis. In 
particular, we set (T ,x p ) — (0.255, 10 -0 51 ) as found by 
Rasia et al. (2004) for their set of simulated clusters. 

In the following, we analyze the implications of the gas 
density and temperature profile of the RTM model on the 
SZE. Contrary to the lensing applications, we limit our 
analysis to the spherically symmetric case, without con- 
sidering any ellipticity in the profiles. The reason for this 
choice is the fact that the gas profiles have been deduced 
by the mean of hydrodynamical (and not simply N - body) 
simulations, so that the hypothesis of axial symmetry for 
the density profile is not enough to assure a similar ellip- 
tical temperature profile. 



The radial dependence of cluster profiles is becoming a 
testing ground for models of structure formation and for 
our understanding of gas dynamics in galaxy clusters. 
Actually, the formation of structures is believed to be 
driven by some hierarchical development, which leads to 
the prediction of self similar scalings between systems of 
different masses and at different epochs. Moreover, the in- 
tracluster gas is generally assumed to be isothermal and 
in hydrostatical equilibrium. From the observational point 
of view, however, the situation is rather controversial and 
yet undeterminated : X - rays observations of poor clusters 
fall belove the self similar expectations, and even if the 
isothermal distribution is often a reasonable approxima- 
tion to the actual observed clusters, some clusters show 
not isothermal distribution ( |Puy et al. 2 000 ) . It turns out 
that the emerging temperature profile is one where the 
temperature increases from the center to some character- 
istic radius, and then decreases again. The central tem- 
perature decrements has been much discussed in terms 
of cooling flows, while the outer temperature decrement 
seems to be confirmed both observationally and numeri- 
cally. On the other hand, the temperature profile described 
in Eq.(@2J reproduces quite well some of these observa- 
tional features: it shows an isothermal core up to 0.2 Rv, 
followed by a steep decrease that reaches a factor two 
lower around the virial radius; the density profiles are self- 
similar roughly s > 0.06, while the gas becomes flatter in 
the inner region. However, these non- trivial clusters pro- 
files need to be observationally tested. The observations 
of the SZE, which are becoming increasingly accurate, can 
be used to probe these properties. Here we analyze the 
radial dependence of the SZE observables, deserving the 
comparison with the observational data in the X - ray and 
SZE domain to a forthcoming paper. 

The temperature shift may be evaluated inserting 
Eqs.lgSJl and into Eqs.O and |J57|| thus obtaining: 



X 1] 



AT = 2k B a T T eo n e 
T m e c 2 
with r\ the so called structure integral defined as : 

r l n e (s) T e (s) 



V 



-dl' 



(46) 



(47) 



which, for a given density and temperature profile, de- 
pends only on the geometry and the extension of the clus- 
ter along the line of sight. In Ea. l|47|) . I is the maximum 
extension of the gas along the line of sight. Measuring the 
lengths in units of the virial radius, it is I = 1 since the 
RTM model is truncated at R v . Note that, usually, one 
takes I — > co thus introducing a systematic bias whose ef- 
fect we will examine later. A simple geometrical argument 
converts the integral in Eq. H47|) in angular form introduc- 
ing the angular diameter distance, Da, to the cluster: 



n = 26 V D A / n e (x)T e ( X )dl' 
Jo 
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Fig. 12. The structure integral r\ vs £ . The upper curve 
has been obtained by fixing I = 1 in Eq. l|47|) . while the 
lower one refers to the case / — > oo . 



26 V D A 



n e (x)T e (x) 



X d X 



(48) 



with x = (I 2 +x 2 + y 2 )/Rl, and £ 2 = (x 2 + y 2 )/R 2 v . The 
integral in Ea. l|48|l can be further simplified introducing 
the auxiliary variable 

. 1 + s 2 



The structure integral n, once evaluated, allows to cal- 
culate the comptonization parameter y, and then the 
temperature shift AT/T, according to the formula l|42l) . 
Adopting a flat ACDM model with (h,ft m ) = (0.7,0.3) 
and typical values for the virial radius and the virial tem- 
perature (R v ~ 2h~ 1 Mpc, T v ~ 8 — 9 keV), we get a 
central number density n e ^ = 5.07 cm -3 and a shift in 
temperature at the cluster centre AT ~ 9.3 mK. 

Let us now investigate in more detail how the peculiar- 
ities of the model affect the structure integral (and hence 
the comptonization parameter and the temperature shift). 
First, we consider the effect of the finite extension of the 
RTM model. To this aim, in Fig. 1121 we plot the structure 
integral r/ as function of £ with the upper (lower) curve ob- 
tained assuming I = 1 (I — > 00) in Ea. H47l) . It is worth not- 
ing that the usual hypothesis of infinite extension may lead 
to significantly underestimate the SZE effect of the model 
by an amount that depends on the value of £. However, 
since what is usually measured is the temperature shift at 
the cluster centre, the relative error is not dramatic being 
less than ~ 10% for £ < 0.2 as it is shown in Fig.^5] Note, 
however, that the error due to the finite cluster extension 
for the RTM model is lower than the corresponding one 
for the standard /3 model ( |Puy et al. 2000| ). 

The temperature profile of the RTM model is approx- 
imately isothermal up to ~ 0.2i?„. Moreover, the isother- 
mality hypothesis is often used in SZE computations. It 
is thus interesting to investigate what is the systematic 
error induced by the simplifying assumption T(r) = T" e .o 
for the RTM model. In this case, the structure integral 
may be analytically expressed in terms of hypergeometric 
functions : 



r) iso (s) = 29 v D A Z(s) 
with 



(49) 



Fig. 13. The relative error e on the structure integral 77 
vs £. It is : e = [r](l -> 00) - n(l = l)]/n(l = 1). 




Fig. 14. The structure integral rj vs £ for the RTM model 
using the correct temperature profile (upper curve) or the 
isothermality hypothesis (lower curve). 
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Fig-E]shows that there is a dramatic change in the struc- 
ture integral (and thus in the SZE temperature shift) 
for the RTM model if we use an isothermal temperature 
profile instead of that given by Ea. H44|) . Actually, even 
if the RTM model has an almost isothermal core up to 
r = 0.2 R v , assuming an isothermal profile leads to a large 
error in the the structure integral. In particular, in the in- 
ner region of the cluster, where the SZE effect is measured, 
rjiso is more than 30% lower than the true rj thus leading 
to a similar error on the temperature shift. Moreover, this 
error is larger than the one due to the finite extension of 
the cluster that dominates in the outer regions. 

8. Comparison with standard gas profiles 

The peculiar non - isothermal temperature profile of 
the RTM model makes it different from most of the 
parametrizations used to describe the gas properties in 
galaxy clusters. It is thus particularly interesting to com- 
pare the SZ signal for the RTM model with that from 
more standard gas profiles. As interesting examples, we 
will consider the j3 and the NFW models. 



8.1. The 13 -model 

The X-ray surface brightness in galaxy clusters 
is commonly fitted using the so called f3- model 
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Fig. 15. The relative error on the SZE signal induced by 
using the best fit [3 - model instead of the correct RTM one 
as function of the dimensionless radius £ = r/R v . 



l|Cavaliere fc Fusco Femiano 1976|) whose density profile 
is : 



- = Pvpx-W ( 1 + ^ 
Pb V X c 



> -M 

2 \ 2 



(50) 



where x c = r c j R v is the core radius in units of the virial 
radius, pop is related to the central electron number den- 
sity n e o = PhPopKurrip) and (3 is a fitting parameter lying 
in the range 0.5 < p < 1. The comptonization parameter 
for such a model is : 



2k B (TT f _ j. k B (TTn e oT e0 
V = 5- / n e T e d/ = 5 77/3 

with 77/3 the structure integral given as : 



(51) 
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(52) 



where we have defined £ = x 2 /x 2 . In order to compare 
the RTM and /? SZE signal we evaluate the quantity 
e y — 1 — yp/vRTM- To this aim, it is convenient to first 
reparametrize the RTM structure integral as follows : 



Vrtm — R v x pi 



1.5 



x+xl X 'o-oo8 (1 + vzr 2,5 ^ 



(53) 



with x Pl = 0.04, x P2 = lO" - 51 and x = x 2 /x 2 pi . With 
these settings, we get : 



= 1- 
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To estimate e y , we have first to choose reasonable val- 
ues for the (3- model parameters. To this aim, we first set 
Pop = po, with p the characteristic density of the RTM 
model, and then choose (P,x c ) by fitting the (3- model 
to the RTM density profile. We obtain as best fit val- 
ues (fl,x c ) — (0.74,0.04) in qualitative agreement with 
(P,x c ) = (0.7, lO- 153 ) given in Rasia et al. (2004). With 
this choice of parameters, the two profiles agree quite well, 
within few %, in the inner regions, while the relative error 
increases up to ~ 10% in the outer zone. Nonetheless, due 
to the radically different temperature profiles, the SZE 
signal is quite different as can be seen in Fig.UHl Actually, 
e y turns out to be significantly different from unity even 
in the region £ < 0.8 where the two models fit each other 
quite accurately. In particular, we find that the SZE signal 
is larger for the RTM than for the /3 - model everywhere 
but in the extreme outer regions of the cluster where the 
situation is reversed. 

Actually, the most interesting implications of the non 
isothermal temperature profile for the RTM model with 
respect the standard (3 -model concerns the detectability 
of a cluster SZE signal, and therefore the statistics of the 
Sunyaev-Zeldovich clusters in different cosmological mod- 
els. Even if all of the physics of the effect is coded in the 
Compton y parameter, it is the total flux density from 
the cluster that is requested from the observational point 
of view. This is found by integrating the comptonization 
parameter over all the cluster face obtaining : 



Y = Jy(6)d 2 6 



(55) 



being 6 the angular position on the sky. Since y is di- 
mensionless, Y is effectively a solid angle. A caveat is 
in order here. Any SZE clusters survey (as for instance 
Plank) has some fixed angular resolution, which will not 
allow to spatially resolve low mass clusters, and even high 
mass cluster can be barely resolved ( Agh anim et al. 1997| ) . 
Therefore, a background y& Compton parameter will be 
present and will be dominated by low mass clusters since 
their higher number density overcompensates their lower 
individual contributions. It is also worth noting that an 
isotropic background would not matter since it could be 
completely removed. Because of this noise, the detectabil- 
ity of a SZ cluster will depend on the average background 
fluctuations that can be estimated as ifBertclma nn 200l)l : 



V 



J y(9')b(9 - 0')d 2 G' 



(56) 



being b(9) the cluster beam profile. A cluster is assumed to 
be detectable if its integrated, beam - convolved Compton 
y parameter is sufficiently large, i.e. : 



Y = / y(6)d 2 e > Y„ 



(57) 



x 2 (l + £l) Jwhere the integral has to be evaluated over the area where 
_ 2 5 \ ~ 1 the integrand sufficiently exceeds the background fluctu- 

dx' \ 



0.13 



(Rations. However, in the following, we will neglect for sim- 
plicity the background noise, i.e. we will simply consider 
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Fig. 16. Dependence of ey on the distance from the clus- 
ter center (in units of R v ). 

Y rather than Y. Even within such a simplified situation, 
we will discuss some aspects of the SZE detectability of 
RTM clusters, which can be easily generalized in presence 
of the background noise. 

If the gas temperature profile is isothermal, the inte- 
grated SZE flux calculated according to Eq.lJSSJ may be 
simply related to the cluster temperature weighted mass 
divided by D 2 A , being Da the angular diameter distance. 
Actually, in an isothermal regime, being dfl = dA/D A , 
Eq. 1|55(1 becomes : 



Y = / y(9)d 2 e cx 



N e (T e ) M(T e ) 



D 2 

A 



D 2 

A 



(58) 



being N e the total number of electrons in the cluster, 
(T e ) the mean electron temperature (which appears in an 
isothermal profile), and M the total mass of the cluster 
(or the gas mass M g = f g M). From Ea. l(5%j) . it turns out 
that an SZE survey detects all clusters above some mass 
threshold which thus has a crucial role for the estimate 
of the SZE cluster counts and its cosmological and astro- 
physical applications fCa rlstrom et al. 2 002). 

However, for the RTM model, the temperature pro- 
file is approximately isothermal only in the inner regions 
(r < 0.2R V ) so that the question of the detectability has 
to be significantly revisited. As a first step, let us compare 
Yrtm, the exact (i.e. computed using the correct tempera- 
ture profile) integrated SZE flux for the RTM model, with 
Yrtm that is evaluated under the isothermal approxima- 
tion. These quantities are given as : 



Yrtm oc 4tt 



A" 



„2.016 



(* pi +z) 2 -5(4 2+a; 4)C,13 



dx 



Y$ M oc 5.33(y - x l J 2 ) 
with : 

+ X [x Pl (2.5x Pl + 1.875A) + 0.375A 2 ] 



y = 



(x pi + x) 2 - 5 



(59) 
(60) 

(61) 



with X = r/R v . Fig.HBlshows how ey = ^ — Y r s ^ m /Yrtm 
increases with the distance from the cluster centre being 
of order 25% at the virial radius. As a result, Eq.JSSJ for 
a RTM cluster is replaced by : 



Y cx 



M eff (T e ) _ vM(T e ) 



n 2 

A 



(62) 
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Fig. 17. Comparison of the integrated SZE flux of a /3 and 
RTM model with ey = 1 — Y^/Yrtm plotted as function 
of the dimensionless distance from the cluster centre. 



being v = Y RTM /Y R S ^ M . Because of the coefficient v > 
1, the SZE detectability of a RTM cluster is naturally 
improved: for a fixed threshold value Yu m , it is possible 
to detect less massive clusters than in the isothermal case. 

As a consequence, a RTM cluster is also detectable at 
lower mass regimes than a (3 - model cluster, as shown from 
a direct comparison. To this aim, we first remind tha the 
integrated SZE flux for a f3 - model is : 



Yp cx 5- x 2 F 1 



2'2 P 



e 



(63) 



In Fig. El we pl°t the relative deviation of the best fit /3 - 
model and the corresponding RTM model. It turns out, 
indeed, that the integrated SZE flux is rather higher for a 
RTM cluster. 

8.2. The NFW model 

Even if it is mostly used to describe the dark matter rather 
than the gas distribution in galaxy clusters, it is nonethe- 
less interesting to compare the SZE predictions for the 
RTM model with the same quantities evaluated for the 
NFW model. 

As a first step, we fit the NFW model to the RTM one 
obtaining cnfw — 5.99 as best fit value for the concen- 
tration parameter. In agreement with Rasia et al. (2004), 
we find that the NFW model slightly overestimates the 
gas density in the range 0.04 < x < 0.4, while it works 
sufficiently well in the outer regions. Assuming an isother- 
mal profile for the gas temperature, the comptonization 
parameter for the NFW model turns out to be : 



Vnfw oc 



(X 2 -x 2 s )(x s + y/l^Tx^)V 



xt - x 2 X 2 



x xUl^xJ^ + xI 1+ VL2+X ' 2 ) x 



In 



(x 2 -* 2 s )yi 



(64) 



with : 



yi = x s (X 2 + x s v^TA 2 ) - Ly/xi - x 2 X 2 , (65) 
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Fig. 18. Relative discrepancy between the Componization 
parameter for a RTM model and its best fitting NFW 
model. 
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having denoted with x s = 1/cnfw> L the cluster length 
along the line of sight and X = r/R v . In Fig.^l we com- 
pare the SZE signal predicted from both models, plotting 
the relative discrepancy e y — 1—Unfw Idrtm between the 
comptonization parameters. We see that the SZE signal 
from a RTM cluster exceeds that from a NFW model as 
yet observed when comparing to the (3 - model. The RTM 
model thus emerges as the most effectively detectable also 
at larger distances from the cluster centre. This is yet more 
clear from Fig.1191 where we compare the integrated SZE 
flux for the RTM and NFW models. 

9. A short comment on possible systematic errors 

There is a potential caveat about the RTM model that 
could affect the main results we have discussed insofar. 
The simulations that have been considered by Rasia et 
al. (2004) to develop the model are based non radiative 
hydrodynamics so that cooling flows and cold blobs that 
may eventually be along the line of sight are not repro- 
duced. Taking into account these effects is a difficult task, 
but we expect, qualitatively, that cooling flows and con- 
duction should lower the central temperature and thus 
increase the SZE signal at the centre. However, we stress 
that detailed simulations also taking into account star for- 
mation and feedback processes are needed to investigate 
how the gas temperature profile (and thus the SZE signal) 
are affected. Some preliminary results have been presented 
UDolag et al. 2004| ), but the effect on the SZE flux has still 
to be investigated. 

As a final remark, it is worth noting that merging of 
clusters has not been considered, but it is likely that this 
does not affect the main results. Actually, the sample of 
simulated clusters analyzed by Rasia et al. (2004) com- 
prises both relaxed, unrelaxed and post - merging systems 
and the RTM model turns out to be a good fit to the 
full sample which is an evidence strongly suggesting that 
merging effects do not alter significantly the cluster struc- 
ture. As a result, the lensing properties of the model are 
likely to be affected only when the merging is in progress in 
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Fig. 19. Relative discrepancy between the integrated SZE 
fluxes for a RTM model and its best fitting NFW model. 

which case an external shear could mimic to first order the 
deviations from spherical symmetry of the outer regions 
of the dark matter halo. A stronger effect is expected for 
the impact of merging on the SZE signal and the X-ray 
emission, since they depend on n eo and n 2 respectively. 
Although further detailed simulations are needed to quan- 
titatively address this question, some partial analytical re- 
sults have been already obtained for some special mergers 
regimes, when the presence of cold fronts marks the late 
merging stages : namely the transonic and the subsonic 
mergers. It turns out that in the transonic regime the fre- 
quency spectrum of the SZE signal g(p) changes, due to a 
shock particle re - acceleration mechanism, depending on 
the concentration, which induces a new electron popula- 
tion. As net effect, the crossover changes still up to ~ 10%. 
In the subsonic case, instead, g(p) remains unchanged, but 
the amplitude of the SZE signal is enhanced in a not neg- 
ligible way, mainly in the interior regions of the cluster 
where it reaches also ~ 30 - 40% IjKoch fc Jetzer 2004jl . 
Moreover, it is well known IjTorri et al. 2004|l that the X- 
ray luminosity overall increases during merging so that it 
is likely that, as net effect, the SZE signal is enhanced 
during cluster merging events. 

10. Conclusions 

Being detectable at high redshift, galaxy clusters are 
promising tools for determining cosmological param- 
eters and testing theories of structures formation. 
Hydrodynamical simulations are able to predict not only 
the dark matter mass distribution, but also the density 
law and the temperature profile of the gas component 
thus allowing a study of both the lensing properties and 
the Sunyaev- Zel'dovich effect due to the cluster. This has 
been the aim of the present paper where we have applied 
this study to the RTM model, recently proposed by Rasia 
et al. (2004) on the basis of the results of a large set of 
high resolution hydrodynamical simulations. 

Assuming spherical symmetry in the density profile, 
we have derived the main lensing properties of the RTM 
model evaluating the deflection angle and the lensing po- 
tential that allows us to write down the lens equations. 
Most of the work has been devoted to a detailed investi- 
gation of the critical curves structure of the RTM model 
since it is the position of radial and tangential arcs (that 
forms just near the critical curves) that gives the most 
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useful constraints on the cluster parameters. The main 
results are summarized below. 

1. The RTM model always gives rise to both radial and 
tangential arcs, but their distance from the cluster cen- 
tre is comparable to that observed in real systems only 
for values of the concentration much higher than typ- 
ically predicted by the numerical simulations of Rasia 
et al. (2004). 

2. Adding a giant elliptical galaxy (described with the 
Hernquist profile) to the cluster potential increases the 
distance between the radial and the tangential critical 
curves and allows to form observable radial arcs for 
lower values of the cluster concentration. 

3. The shape of the critical curves may be significantly 
altered by a perturbing shear (mimicking an inter- 
nal ellipticity or the effect of tidal perturbations) thus 
changing the positions of both radial and tangential 
arcs by an amount that can be up to <~ 30% for a 
shear strength 7 = 0.15 (see Table 1). 

4. Fitting an NFW model to a cluster that is intrinsically 
described by the RTM model leads to overestimate the 
cluster mass by an amount that depends on the con- 
centrations crtm and cmfw used in the modeling. 

The work presented should be complemented by a detailed 
study also taking into account the presence of substruc- 
tures (predicted by the CDM paradigm of structures for- 
mation) in the dark halo cluster. The total lensing poten- 
tial thus should be made of the sum of the contributions 
from the elliptical RTM model, the external shear (due, 
e.g., to tidal perturbations or large scale structure), the 
bright cluster galaxy and the satellite haloes. This ap- 
proach is quite complicated given the high number of un- 
known parameters entering the modeling, but one could 
resort to X - ray data to constrain (at least) the RTM pa- 
rameters (M v ,x p ). However, the most compelling test is 
a direct comparison with real systems showing tangential 
and radial arcs. To this aim, adding the BCG contribution 
(and eventually the shear term) to the RTM lensing po- 
tential should give a sufficiently accurate cluster model to 
be compared with the data on the arcs position. Having 
determined the galaxy scalelength and surface brightness 
from photometric observations, this method should allow 
us to test whether the RTM model may reproduce the ob- 
served arcs positions and to constrain its parameters. This 
will be performed in a future paper. 

One of the most interesting feature of the RTM model 
is that the temperature profile is approximately isother- 
mal only up to 0.2 R v . We have investigated the impli- 
cations of the RTM model for the intracluster gas on the 
SZE evaluating the structure integral that determines the 
temperature decrement. We have thus taken into account 
both the finite extension of the model and its peculiar 
temperature profile estimating the errors induced by the 
usually adopted simplifications of infinite extension and 
isothermal temperature. The main results are as follows. 



1. Neglecting the finite extension of the cluster systemat- 
ically underestimates AT/T by an amount that is less 
than 10% in the inner regions of the cluster so that it 
may be neglected in a first order analysis. 

2. Using T = T e fl as temperature profile instead of that 
found by RTM leads to underestimate AT/T up to 
30% (in absolute value) in the centre. 

3. The comptonization parameter y for the RTM model 
is higher than that of both the (3 and the NFW mod- 
els even if the parameters are chosen in such a way 
that the gas density is well fitted by the three mod- 
els. Using the best fitting (3 (NFW) model instead of 
the correct RTM one underestimates y up to <~ 40% 
(~ 52%, respectively) in the inner cluster regions. 

4. The non isothermal temperature profile leads to an 
integrated SZ flux which is higher for the RTM model 
than for both the (3 and NFW models by an amount 
that depends on the distance from the cluster centre, 
but can be as high as ~ 45% and ~ 90% for the (3 
and NFW models respectively. As a result, less massive 
clusters should be detected in SZE survey if the RTM 
model is indeed the correct one. 

As a final remark, we would like to stress that, in our 
opinion, a combined analysis (from the theoretical and 
observational point of view) of the both the lensing prop- 
erties and the SZ temperature decrement could be the best 
method to validate a given cluster model independently on 
the peculiarities of the numerical simulations inspiring it. 
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